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ABSTRACT 

We present the results from a survey of NGC 2244 from 3.6 to 24 fim with the Spitzer Space 
Telescope. The 24/j,m-8/im-3.6^m color composite image of the region shows that the central cavity 
surrounding the multiple O and B stars of NGC2244 contains a large amount of cool dust visible 
only at 24/im. Our survey gives a detailed look at disk survivability within the hot-star-dominated 
environment in this cavity. Using mid infrared two color diagrams ([3.6]-[4.5] vs [5.8]-[8.0]) we iden- 
tified 337 class II and 25 class I objects out of 1084 objects detected in all four of these bands with 
photometric uncertainty better than 10%. Including the 24 /Ltm data, we found 213 class II and 20 
class I sources out of 279 stars detected also at this latter band. The center of the class II density 
contours is in very good agreement with the center of the cluster detected in the 2MASS images. We 
studied the distribution of the class II sources relative to the O stars and found that the effect of high 
mass stars on the circumstellar disks is significant only in their immediate vicinity. 
Subject headings: open clusters: individual (NGC 2244); stars: pre-main sequence; (stars:) planetary 
systems: protoplanetary disks; infrared: stars 



1. INTRODUCTION 

The effect on star formation of an environment domi- 
nated by high mass stars is still an unresolved problem 
of astrophysics. Strong stellar winds from high mass O- 
, and B-type stars can trigger star formation by com- 
pressing and the interstellar materia l, which eventually 
becom es gravitationally unstable (El megreen fc Ladal 
Il977f h Also photoevaporation can heat the surface of 
an interstellar cloud and lead to a radiatively driven im- 
plosion, thus playing a similarly (if not more) impor- 
tant role in triggerin g star formation in this environment 
(| Adams et al.l 120041) . On the other hand the strong ex- 
treme ultraviolet (EUV) and far ultraviolet (FUV) radi- 
ation from the same hot stars might so vigorously pho- 
toevaporate the material around the forming stars that 
it truncates the star formation process. A better under- 
standing of the photoevaporation of circumstellar mate- 
rial would therefore be an important step toward clari- 
fying the overall effects of hot stars upon the process of 
star formation. 

The discovery of " proplyds" in Orion and other 
high-mass sta r -form ing regions (e.g. lO'dell et al.l Il99l 
lO'dell fc Wenl I1994D led to a number of theoretical 
studies of the photoevaporation of young circumstel- 
lar disks by external radiation fields from neighboring 
hot sta rs (Ijohnstone et al.lll998t IRichling fc Yorkdll998l 
20001: IMatsuvama et all l2003t [Hollenbach fc Adams! 
2004t iThroop fc BalM 120051 ). These studies have shown 
that the outer parts of disks evaporate quickly, produc- 
ing structures that match the appearance of the pro- 
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plyds well. iJohnstone et al.l (|1998l ) presented a model 
for the photoevaporation of circumstellar disks by an ex- 
ternal source of UV radiation. They applied this model 
to 8 1 Ori C and concluded that circumstellar disks are 
rapidly destroyed by the external UV radiation field. Ex- 
cluding viscosity effects, they predict final disk sizes of 
r < 1 AU, at distances within 0.3 pc of the illuminating 
source and reached about 1 million years after the onset 
of the UV ill umination of the disk. A number of other 
authors (e.g. IRichling fc Yorkd 120001: IMatsuvama et al.l 
l27)03llHollenbach fc Adamsll2004l : iThroon fc Ballvll2005h 
predict similar (or shorter) timescales but differ regard- 
ing the extent of the disk erosion. 

InfraRed Array Camera (IRAC) and Multiband Imag- 
ing Photometer for Spitzer (MIPS) observations of high 
mass star forming regions are perfect tools for testing 
these theories. They probe the inner disk emission out 
to «10 AU and can cover large areas containing enough 
stars to allow the use of statistical methods. Here we 
present results from Spitzer imaging of a massive young 
cluster, NGC 2244 

NGC 2244 is a young, nearby open cluster contain- 
ing 7 O and 65 B stars. These high mass members of 
the emerging cluster are probably responsible for the 
evacuation of the central part of the spectacular HII re- 
gion, the Rosette Nebula. The cluster has been studied 
at optical and X-ray wavelengths (e.g. lOgura fc Ishid"al 
198ll:lLi et a"iT2 002: Li 2003: iGregorio-Hetemet 1 _aL|1998l : 
Berghofer fc Christianll2002l : iPerez et allll987l 119891 ). It 
is among the most suitable areas for studying active star 
formation in a high-mass environment. 
griekeSNsGiyzdi^^d^oabgtts.ArfeAiJ.ckpc from the Sun and 
has an age of about 2-3 Myr dPark fc Sund [2 002: 
Ogura fc Ishida|[l98l IPerez et al.lll987tlHensberge et all 
20001 : IPerez! Il991h . Most of the age determinations 
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Fig. 1 .— False color image of NGC 2244. B:3.6 /he, G:4.5 (an, R:8.0 fira. The image is centered at RA=06 h 32 m 18 s and DEC=04°51'46' 
covering about 30'x30' area. 



used color-magnitude diagrams in differe nt filter-systems 
and w ith slightly different approaches. lOgura fc Ishida 
(1981) used the turn off point of the dereddened UBV 
color-magnitude digram and the posit ion of the ea rliest 
O stars to estimate an a ge of 4±1 Myr. [Peredd 19911) used 
the statistical method of lSchroeder &: Cominsl (|1988t ) and 
derived 2.1 ± 1.6 Myr. Probably the most accurate age 
and di stance estimate was performed bv lHensberge et all 
(2000). They used uvby photometric and high resolu- 
tion spectroscopic observations of the cluster member 
and eclipsing binary V578 Mon. They applied a Fourier 
disentangling technique resulting in the distance and age 
1. 39 ±0.1 kpc and 2.3 ± 0.2 Myr respectively. 

iPark fc Sungj (|2002f) calculated the ages of the pre- 
main sequence (PMS) and main-sequence cluster mem- 
bers (identified from the combination of photometric, 
spectroscopic, and p roper motion studies) using the evo- 
lutionary models of Schall er~et alj (I1992D for mai n se- 
quence stars and thelD'Antona fc MazziteTil (|1994h and 
ISwenson etHI (|1994h PMS models. The turnoff age of 
the cluster was 1.9 Myr while the mean age of the low 
mass PMS stars (2.5 M > M > 0.5 M Q ) was under 1 
Myr with an age spread of 6 Myr, although the small 
numbers of PMS stars and low photometric quality for 
these objects makes the estimate of the age spread rather 
unreliable. 

In this paper we present Spitzer IRAC and MIPS 24/im 
observations centered on the main OB core of NGC2244. 
We describe our observations and data reduction tech- 
niques in §2. Our results and analysis are presented in 



§3 and §4, followed by our conclusions in §5. We base 
our analysis on an assumed age of 3 Myr and distance of 
1.5kpc. 

2. OBSERVATIONS 

Observations of N GC2244 were ob tained on 2004 
March 09 with IRAC (|Fazio et al.ll200l . The IRAC sur- 
vey covers about half a square degree centered on the 
main OB core. The 12s high-dynamic-range mode was 
used to obtain two frames in each position, one with 

0. 4s exposure time and one with 10.4s. The observa- 
tion of each field was repeated twice with a small off- 
set that allowed 20.8s integration time for each pixel. 
The frames were processed using the SSC IRAC Pipeline 
vl3.2, and mosaics were created from the basic calibrated 
data (BCD) frames using a custom IDL program. 

Source finding and aperture photometry were car- 
ried out using Phot Vis version 1.10 which is an 
IDL- GUI based photo metry visualization tool. See 
iGutermuth et alj (|2004D for further details on Phot Vis. 
The radii of the source aperture, and of the inner and 
outer boundaries of the sky annulus were 2.4, 2.4 and 
7.2 arc-second respectively. The calibration was done by 
using large aperture measurements of standard star ob- 
servations. The zero point magnitudes of the calibration 
were 19.6600, 18.9440, 16.8810, and 17.3940 for channel 

1, 2, 3, and 4 respectively. Aperture corrections of 0.21, 
0.23, 0.35 and 0.5 mag were applied for channels 1, 2, 3, 
and 4 to account for the differences between the aperture 
sizes used for the standard stars and for the NGC 2244 
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Fig. 2. — False color image of NGC 2244. B:3.6 fim, G:8.0 (im, R:24 fim. The center and the covered area are the same as in Fig[T] 



photometry. 

The detected sources were examined visually in each 
channel to clean the sample of non-stellar objects and 
false detections around bright stars. We accepted as good 
detections those with photometric uncertainties less then 
0.1, which allowed limiting magnitudes of 17.0, 16.4, 14.2 
and 13.6 at 3.6 /im, 4.5 /im, 5.8 /im, 8.0 /im respectively. 
We detected more than 10000 sources at 3.6 /im, 8000 
at 4.5 /im, and over 2000 and 1000 in the 5.8 /im and 
8.0 /mi images respectively. There are 1084 sources that 
were detected in all 4 channels, a criterion that elimi- 
nated almost all of the stars fainter than 15.0 mag at 
3.6/im. We rejected the remaining eight sources fainter 
than 15.0 mag to help minimize the contamination from 
background galaxies. Based on the 2MASS magnitudes 
of the detected sources, the limiting mass of our survey 
is around 0.8- 0.9 M w . 

The MIPS (jRieke et al.l I2004H 24 /im survey was ob- 
tained on 2005 May 16. It covers about 0.8 square de- 
grees, including the region of the IRAC survey. The 
observations were carried out in scan map mode. The 
frame s were processed u sing the MIPS Data Analysis 
Tool (jGordon et al.ll2005l ). The high 24 /im background 
towards NGC2244 made automatic source extraction ex- 
tremely difficult so we visually identified the sources on 
the background subtracted images. We identified 501 
24 /im sources; most of them are concentrated in the 
area where the O stars are located. PSF fitting in the 
IRAF/DAOPHOT package was used to obtain photome- 
try. Our limiting magnitude at 24 /im is 10.7 using a 7.3 
Jy zero-point for the 24 /im magnitude scale. 279 of the 



24 /xm sources were also detected in the IRAC frames 
with acceptable uncertainty. Almost 45% of the 24/im 
sources have no counterparts within the IRAC images. 
There are two reasons: 1.) the 24/im survey covers a 
larger area than the IRAC survey (causing loss of more 
than 150 sources); or 2.) the photometric uncertainties 
of the 24/im sources are sometimes too high in the IRAC 
images so they are eliminated by the error cut. It is also 
possible that some of the 24/im sources with no counter- 
parts are spurious detections due to the complex back- 
ground, but careful visual examination of the sources lim- 
ited this possibility only to a couple of cases. Many of 
these detections are probably high-redshift galaxies. 

Our final photometric catalogue with sources detected 
in all four IRAC channel sorted by right ascension is 
presented in Table 1. 

3. RESULTS 

3.1. Extended structure 

False color images of the NGC2244 region are shown 
in Figure[T](B:3.6 /im, G:4.5 /im, R:8.0 /im) and[2](B:3.6 
/im, G:8.0 /zm, R:24.0 /im). The prominent central cav- 
ity of NGC2244 can be identified in both images, though 
it is partially filled with dust emitting at 24 /im. The 24 
/xm emission is extremely strong in an elongated struc- 
ture in the middle of the main OB core. One interest- 
ing feature of this elongated structure is that it has a 
bowshock-like geometry facing the nearest O and B stars. 
The highest mass star in the region is HD46223 (spectral 
type 05), which is situated southwest from the brightest 
star, HD46241 (K0V foreground object). There is a bow 
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Fig. 3. — Color images of the vicinity of HD46150 8fim. Left panel:R,ed 24^tm; Green Blue 3.6/jm. Right panel: Red Green 

4.5/zm; Blue 3.6/Ltm. The centers of the images are at RA=06 h 31 m 55 s DEC=04°56'22", the covered area is 75"x62". 



TABLE 1 

The first five rows of the final photometric catalogue of sources detected in all four IRAC channels. The whole 

CATALOGUE IN ELECTRONIC FORMAT IS AVALILABLE AT THE ApJ WEBSITE. 



RA(J200) 


DEC(J2000) 


[3.6] 


a([3.6]) 


[4.5] 


,7(4.5]) 


[5.8] 


<r([5.8]) 


[8.0] 


<r([8.0]) 


[24] <x([24]) 


[dcg] 


[dog] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] [mag] 


97.8174 


5.0398 


14.254 


0.019 


13.844 


0.020 


13.548 


0.098 


12.675 


0.092 




97.8195 


5.0555 


12.380 


0.008 


12.358 


0.009 


12.298 


0.031 


12.298 


0.050 




97.8256 


5.0197 


12.311 


0.006 


12.285 


0.008 


12.320 


0.027 


12.157 


0.053 




97.8291 


4.9182 


12.075 


0.006 


12.037 


0.008 


11.976 


0.027 


12.069 


0.082 




97.8298 


5.0017 


11.605 


0.004 


11.461 


0.005 


11.450 


0.016 


11.332 


0.038 





shock-like 24 /im feature associated with HD46223 but it 
does not seem to be as strong as the one associated with 
the other high mass 06 star, HD46150, and several B 
stars just south of the elongated 24 /jm structure. These 
structures may be a result of stellar winds interacting 
with the dust. 

Another interesting feature on Figure [2] is a strong 
24 Lim point source with an extended tail right next to 
HD46150 (see F igure El for the en larged images). In a 
separate paper (jBalog et al.ll2006T ). we show that it is 
probably a class II source with a photoevaporating disk. 

There are two other remarkable objects in the images. 
One of these is a small patch on the 8 /xm image extend- 
ing north from a class II source (see Figure [4] left panel 
for the enlarged image and Section 3.2 for the description 
of the classification of the sources). It seems to harbor 
a point source visible only at 8 ^m. The structure looks 
like a small globule in the optical images. Very weak 
extended 24/im emission can be detected in association 
with this 8 /jm feature. The second is a formation similar 
to a "pillar of creation" to the southeast from the main 
24/zm structure. It has a class II object at the tip, and 
there are a class I and several class II sources to the west 
of this formation (Figure 2] right panel). This object has 
a linear size of 2.2 pc which is very similar to the large st 
of the "pillars of creation" in M16 ([Hester et al.lll996f) . 

3.2. Young Stellar Object Classification 

The young stellar objects (YSOs) can be charac- 
terized by their spectral ener gy distributions (SEDs) 
(jAdams et al.lll987t lLadalll987t ). Based on the SED we 



can distinguish 4 main types of YSOs: class 0, I, II, 
III. The class 0/1 objects are deeply embedded and their 
SEDs peak in the submillimeter or the far infrared indi- 
cating that the source of the emission is cold dust. The 
class II sources are optically visible stars with infrared 
excess emission that is attributed to a disk surrounding 
the central object. Class III sources have almost no in- 
frared excess and their photometric properties are very 
similar to normal main-sequence stars. We can use a 
mid-infrared two-color diagram to distinguish between 
class I, II and III objects. 

We identify 337 class II and 25 class I sources in the 
IRAC two color ([3.6H4. 5] vsf5.8H8 . 0l) di agram using 
the method described bv lAllen etall (|2004h (Figure [5]). 
There are also 14 sources with uncertain evolutionary 
status Q3.6]-[4.5] < 0.4 and [5.8]-[8.0] > 1.1). For them, 
the large value of the [5.8]-[8.0] color suggests class I 
status while the small [3. 6] -[4. 5] color would indicate 
a class II source. We adopt the classification scheme 
of iMegeath et al.1 (|2004h and assign class I/II status to 
these stars. We double checked our classifications using 
other color-color and color-magnitude diagrams like J- 
H vs K-[24] or J-H vs H-K from the 2MASS data. In 
each diagram the class II and class I sources occupy the 
regimes predicted by model calculations or by observa- 
tions of other class II/I sources. 

Another method of identifying stars with infrared ex- 
cess (class I or class II sources) is to combine MIPS and 
IRAC photometry (Figure [6]) . This method should be 
more robust than the previous one because reddening 
hardly affects the 24 /im observations. Also, the separa- 
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Fig. 4. — Color images of the 8^tm patch (left panel, red: 24/xm; green: 8/xm; blue: DSS Red, the center of the image is at RA=06 h 31 m 43 m 
DEC=05°03'10", the covered area is 58"x67") and the "pillar of creation" (right panel, red: 24/xm, green: 8^tm, blue: 3.6/mi,the center of 
the image is at RA=06 h 33 m 10 s DEC=04°45'55", the covered area is 252"x296". 
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Fig. 5. — The IRAC color-color diagram. Dots: "normal" 
stars; filled circles: class II sources; open triangles: class I sources; 
open diamon ds: class I/II sourc es. The classification criteria were 
adopted from lAllen et al.l H2004T ) 

tion of stars with different evolutionary status is larger 
on the combined IRAC-MIPS two color diagram. How- 
ever the 24/im observations are less sensitive and so not 
as complete as the IRAC-only sample. Using the [3.6]- 
[4.5] vs [8]-[24] color-color diagram and the color criteria 
defined by Muzerolle et al. (in preparation) we found 20 
class I and 213 class II sources out of 279 24/mi sources 
detected in the region covered by both IRAC and MIPS 
observations. 

We note that there are some cases where sources identi- 
fied as class I on the [3.6]-[4.5] vs [8.0]-[24] plots are clas- 
sified as class II sources on the IRAC two color diagram 
and vice versa. This indicates that the adopted boundary 
between class I and class II sources is uncertain and in- 
troduces some ambiguity in distinguishing the two types 
definitively based on their infrared colors. Where the 
source showed class II nature on the [3. 6] -[4. 5] vs [8.0]- 
[24] diagram but class I or class I/II on the IRAC color- 
color diagram, we classified it as class II because these 
objects are probably highly reddened class II sources. 
Adding these sources to the IRAC class II sample would 
modify our results to 362 class II sources and 12 class I 



1.5 










1 


MS 


class II 




class 1 


0.5 




■•■ v 


J . 

































0.5 












1 1 


2 3 


4 


5 



[8.0R24] 

Fig. 6. — Combined IRAC-MIPS color-color diagram. Dots: 
"normal" stars; filled circles: class II sources; open triangles: class 
I sources; open diamonds: outliers (stars with uncertain classifica- 
tion). The classification criteria were adopted from Muzerolle et 
al. (in preparation). Solid lines show the boundaries of each class 



sources. 

A small number of young stars might be misidentified 
background planetary nebulae or AGB stars or galax- 
ies. The effects of thes e contaminations are usually small 
(jMegeath et al.ll2004f ). The method of Gutermuth et 
al. (in preparation) uses different color-color and color- 
magnitude diagrams to estimate the number of contami- 
nating sources. It shows that the number of sources that 
might be AGN or PAH emitting galaxies is below 1%. It 
is also possible that a 24/im source invisible at shorter 
wavelengths is projected onto a pure IRAC class II source 
confusing the classification methods. This configuration 
might appear as a class I object showing a large excess 
at 24^im on the [3.6]-[4.5] vs [8.0]-[24] plot while on the 
IRAC two color diagram it seems to be a normal class 
II object. The probability of such an event is extremely 
low (P < 10 -6 ). Chance coincidences with background 
24 /^m-emitting galaxies are also possible, but unlikely in 
the relevant flux density rang e for class I sources, g iven 
typical galaxy number counts (jPapovich et al.|[2004f ) and 
our flux limits. 
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Fig. 7. — 24 (im image with stars with IRAC color excess (class I sources: red; class II sources: green; class I-II: yellow) and O stars: 
blue overlaid 



4. DISCUSSION: THE DISTRIBUTION OF IR EXCESS 
SOURCES AND ITS IMPLICATIONS 

We show the 24 fj,m image with the class II and I 
sources and the O stars overlaid in Fig The class 
II sources are concentrated mainly in the OB core. We 
used smoothed stellar density contours to determine the 
center of the class II population (Figure [5]). Our calcula- 
tion shows that the center of the concentration of class II 
sources is at RA=6 h 31 m 58 s DEC=4°54'51" (just south- 
west from the main 24 ^m emission peak), which is in 
very good agreement with the coordinate of t he cen ter of 
the cluster derived from 2MASS data byO (|2005l) . 

4.1. Circumstellar disks 

The left panel of Figure [5] shows the number of stars 
vs their projected distance to the closest O star in equal 
area (0.8 pc 2 ) bins. The hollow histogram represents 
stars without excess while the filled histogram represents 
class II sources. The right panel of the same figure shows 
the fraction of stars with excess vs their distance to the 
closest O star. Given that the infrared excess is gener- 
ally attributed to radiation from a circumstellar disk, we 
consider each excess source a star with disk. There is 
an excess of sources without disks in the first distance 
bin (42 without disk vs 16 with disk) very close (d < 
0.5 pc) to the nearest O star. In the other bins the disk 
ratio seems to maintain a constant level. The number of 
non-disk stars as a function of projected distance from 
the nearest O star as shown in the left panel of Figure 9, 
levels off around 2.5 - 3.0 parsec. This suggests that the 
stars outside this radius are probably non-members and 
belong to the foreground or background population. 

We can divide our data into three subsamples: 1.) the 



first bin very close to the O stars, 2.) between 0.5 and 2.5 
pc from the O stars where the population of the cluster is 
significant; and 3.) farther than 2.5 pc from the O stars 
where the background and foreground population takes 
over. We subtracted the average background population 
(1 source/bin in the disk and 3 sources/bin in the non- 
disk population) from each bin and calculated the disk 
ratio in the two other regions. The disk ratio is about 
27% if the distance from an O star is less than 0.5 pc. In 
the remaining region where the cluster population is sig- 
nificant the disk ratio is about 45.4% ± 6%. The overall 
disk fraction in the cluster is 44.5%. This is somewhat 
smaller but in reas onable agreement with the results of 
lHaisch et al.l (|2001h . who found a typical disk fraction 
~50% for 3 Myr-old region s using K-L exce s s as a n indi- 
cator of circumstellar disks. iMamaiek et al.l (|2004f ) found 
the same disk fraction at 3 million years using N-band 
excesses. The bin to bin fluctuation is on the order of 
5-10% in agreement with the average statistical error of 
the disk ratio. The small variation with radius indicates 
that the disk fraction is not correlated with the distance 
from the O stars outside a certain distance. However, 
close to the O stars we observe a drop of the disk frac- 
tion. To calculate the significance of this drop, we used 
the binomial theorem. We assumed that the sources close 
to and far away from the O-stars were drawn from the 
same parental distribution, that is, that the average pro- 
portion of disks was the same in both samples and that 
the apparent difference was due to statistical fluctuations 
distributed as in the binomial distribution. We adjusted 
the common value of the average to maximize the proba- 
bility, and found that there was less than a 2the O stars 
having the same intrinsic average proportion of disks as 
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offset RA [arcmin] 

Fig. 8. — Smoothed stellar density contours of the cla ss II sources of NGC2244. The (0;0) coordinate corresponds to the center of the 
IR cluster detected in the 2MASS database by Li (2005). red stars denote the known O stars detected in our survey. 



for the far region. This behavior may be consis tent with 
the model calculations of lJohnstone et all (|1998f) showing 
that a disk can be destroyed in a couple of million years if 
it gets too close (< 0.3pc) to an 06 star. However, there 
is only a trend against disks near the O stars, not a com- 
plete lack of them. Four out of five O stars in our field 
have at least one confirmed class II source within this 
radius, including the two highest mass O stars, HD46223 
(05) and HD46150 (06). Altogether we see eight class 
II sources within a projected radius of 0.3 pc from one 
of the five O stars. Including the photoevaporating disk 
candidate mentioned in section 3.1, we see a net of nine 
class II sources in the close vicinity of an O star. 

We examined the possibility that these sources just 
seem to be close to the O stars and actually are in front of 
or behind them outside the 0.3 pc radius. We performed 
a simple Monte Carlo simulation placing 5 O stars and 
362 class II sources in a sphere with a radius of 6.5 pc 
and counted the class II stars that were inside the 0.3 
pc projected radius of any O star. After repeating the 
experiment 100000 times, to get accurate statistics, we 
found the most probable number of false companions is 
3 or less (more than 50% of all the cases). The proba- 
bility of seeing 8 or more fake companions is under 3%. 
We repeated the calculations using a 0.5 pc radius to 
reduce the effect of low number statistics and to com- 
pare our calculation with the results shown on Figure [9l 
In NGC 2244 we see 16 class II sources within a 0.5 pc 
radius of an O star. In this case our experiment shows 
that the most probable number of chance alignments is 9. 



The probability of seeing 16 or more stars in this region 
is about 4%. 

Given the results of the simulations we can rule out 
the possibility that all of our class II sources are out- 
side the radius of influence of the O star and are only 
chance alignments. To explain the presence of class 
II sources we consider the scenario that the low mass 
star is just passing by the O star and recently en- 
tered the "danger zone" , making the infrared excess a 
short lived phenomenon. In this case, given t he average 
veloci ty dispersion in the clu s ter (^ 10km/s, iLiu et ali 
(1989) or 6.8 km/s. ILiu et~aTl (|199lD ) and assuming ra- 
dial orbits, a low mass star would spend approximately 
~ 3 — 4 x 10 4 yr in t he vicinity of th e O star as it 
passes. (The sample of ILiu et al.1 (119891) con tains only 
3 stars while the sample of ILiu et al.l (|1991l ) contains 
only 8 stars with sensible radial velocities, which makes 
the ra dial velocity dispe r sion rather uncertain.) Mod- 
els by iRichling fc Yorkd (|1998l I2000D predict that this 
circumstance will expose a large amount of dust that 
is entrained in the evaporating gas, leading to bright 
IR emission. The Orion proplvds (lO'dell et al.l 119931 : 
lO'dell & We"n1ll994HSmith et al.ll2005[) are examples, and 
Figure [3] shows a similar object in NGC 2244. They can 
be identified in MIPS images by long "cometary" tails of 
warm dust as the photon pressure ejects the grains from 
the vicinity of the hot stars. 

Although the majority of the class II sources in the 
vicinity of the O and B stars do not have the expected 
morphology they might still be "cometary" sources seen 
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Fig. 9. — The number of excess (filled histogram) and non excess (hollow histogram) sources (left panel) and the disk ratio (right panel) 
vs projected distance from the nearest O star. The solid horizontal line in the right panel shows the average disk ratio of 45.4% between 
0.5 pc and 2.5 pc. The values of more than 100% arise when the number of non-excess sources is below the background level in some of 

Vn 

N 



the bins. The errorbars in the right panel represent 1 a values calculated through a 



in each bin. 



head on or tail on. The modeling of the object in Figure 
El shows that its t ail lies nearly on the plane of the sky 
jBalog et al.H2006l) . At the distance of NGC 2244 about 
3-5 out of 6 such sources would be unresolved with MIPS 
given the 6" beam size and assuming a random distribu- 
tion of tail orientations. 

To sum it up, chance alignments (~9±1) and short 
lived unresolved cometary structures (~11±3) would ac- 
count for most of the class II sources very close to the 
O stars (d < 0.5). The drop in the portion of class II 
sources also suggests the presence of some objects with 
truncated disks. 

Outside a certain radius (0.5 pc in our case) the disk 
fraction does not seem to depend on the distance from 
the O stars. The change in disk properties near 0.5 pc 
may let us test models of photoevaporation. For this 
regime we can consider the limiting case of the disks be- 
ing long-lived. In the current photoevaporating mod- 
els the behavior of the inner zones (within 50 AU of 
the central star) differs substant i ally a mong the stud- 
ies. For example. Uohnstone et al.l (|1998f ) calculated that 
the disks evaporate quickly down to a radius of about 20 
AU, and thereafter more slowly to one of about 2 AU. 
They argue that effects of viscosity not included in the 
calculations would probably result in complete destruc- 
tion o f the disks. The calculations of Richli ng fc Yorkel 
(2000) generally show the evaporation slowing substan- 
tially in the 50 to 100 AU region (but they only follow 
the evolution to 10 years). iHollenbach fc Adamsl (12001 ) 
and iThroop fc Ballvi (|2005l ) calculate that the evapora- 
tion slows substantially at a radius of 5 to 10 AU, requir- 
ing 10 5 years in the lat ter case and slightly long er in the 
former. The models of Matsuvama et all (|2003 ) require 
a few times 10 5 years to evaporate to a limiting radius 
of about 2 AU. 

Some studies (see e.g. Uohnstone et al.l 119981 : 
iRichling fc Yorkel I1998D derive a power law relation 
between the photoevaporation rate and the distance 
from the O star (M p h ~ d~ a ). The exponent of this 
function (a) is in the range of about 0.7 and 1.2. These 
quantitative predictions are amenable to test with our 
observations. To estimate the effect of photoevaporation 
on the disk fraction and to see which slope is most 
probable in the case of NGC2244, we performed a Monte 
Carlo simulation. 



Similar to the previous case where we examined the 
probability of chance alignments, we scattered 5 O stars 
in the volume of NGC2244. Then we let a star with 
the disk randomly walk among the O stars with a veloc- 
ity of 10 km/s. We randomly assigned an initial disk 
mass between 0.001 and 0.1 Mp, to the s ource based 
on the study of I Andrews fc WiUiamsl (|2005h . who found 
that the distribution of disk masses in a similar mass 
range probed by our observations is roughly equal and 
independent of stellar mass. We calculated the pho- 
toevaporation rate in every 0.05 pc step based on the 
distance from the nearest star and the photoeyapora - 
tion rate - distance relation of Rich ling fc Yorkel (|1998f ) 
(M ph = 1.46 x 10- 6 Af o yr- 1 ( TB ^U r )- Q ). They consider 
cases that yield a=0.767 or a=1.134; we therefore have 
let a vary between these values. The output of our sim- 
ulation was the total mass loss and the difference be- 
tween the lost mass and the initial disk mass after 3 
million years. If the latter was less than zero we con- 
sidered the disk totally evaporated. This very simple 
model does not take into account that the photoevap- 
oration may slow down after reaching a certain r adius 
(|Hollenbach fc Adamsl [200l IThroop fc Ballvi I2005D and 
also ignores the effect of gravity (when we calculate the 
trajectory of the stars), the cumulative effect of more 
than one O star's radiation field and the central star's 
own UV radiation field. Other than the first, the ig- 
nored effects would speed up the photoevaporation, so 
our calculated disk survival rates are upper limits (un- 
less in fact photoevaporation slows at some radius > a 
few AU). 

To get reasonable statistics we ran the simulation 
10000 times for the same O star configuration chang- 
ing the initial mass of the disk and the initial position of 
the disk-bearing star. Then we repeated the experiment 
several times for different O star configurations. We also 
considered two extreme cases: 1.) the disk-bearing stars 
could freely leave the cluster; and 2.) the disk-bearing 
stars bounce back from the cluster boundaries. We found 
that the difference between the two cases is within the 
error of our simulations, which is on the order of 2-3%. 

Giv en th e fact that the ve locity dispersions of lLiu et al.l 
(fl989h and lLiu et all (fl99lh are very uncertain and much 
larger than the virial velocity (~ 1 km/s) of the cluster, 
we repeated the experiment using the velocity disper- 
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sion of 2 km/s measured in the Orion Nebula Cluster 
(jSicilia-Aguilar et al.ll2005f ). The effect of this velocity 
change on the disk survival (~1% for each slope we stud- 
ied) is smaller than our errors so we use our original 
values in the following. 

Assuming that a 3 million year old cl uster without 
O stars has a disk fractio n around 50% (|Haisch et al.1 
I200H iMamaiek et ~aTl l2004D. we found th at the value of 
a = 0.767 suggested by Richl ing fc York! £l998) resulted 
in too few surviving disk (5%). As we make the slope 
steeper we see that the disk fraction initially rapidly rises 
(36% with a = 1.05), and then the rise slows down giv- 
ing a disk fraction of 41% at a — 1.134. In spite of the 
weaknesses of our simple model, we can clearly conclude 
that the data favor the steeper slope (a > 1) in the pho- 
toevaporation rate vs distance relation. Together with 
the small difference (< 10%) between the disk fraction 
of clusters with and without O stars and the fact that the 
disk fraction is constant outside a certain radius from the 
O stars, our lower limit a > 1 indicates that the effect 
of high mass stars is significant only in their immediate 
vicinity (d < 0.5 pc). 

4.2. Class I objects 

Many of the class I sources are located around the 
rims of the surrounding nebula, which might indicate 
self-propagating star formation triggered by the expand- 
ing shell driven by OB stars in the core of the cluster. 
However some class I sources appear to be in the evacu- 
ated region around the O stars. It is possible that some 
of these sources are only projected onto the cavity and 
actually they are located behind it. If we consider a con- 
figuration in which the shape of the interface between 
the denser gas and the evacuated region is a half-sphere 
behind the cluster and the sources are evenly distributed 
on this interface, we find that projection can account for 
all of the sources apparently in the evacuated region. 

The fact that probably all class I sources are close to 
their birthplaces, the interface between NGC2244 and 
the surrounding gas, suggests that these objects are re- 
ally much younger than the class II sources that are 
spread throughout the whole cluster. On the other hand 
class I sources represent a group of stars formed by trig- 
gered star formation so there might be young class II 
sources associated with them. Further observations are 
needed to study this problem in more detail. 

The ratio of class I/class II sources is about 7% in the 
IRAC-MIPS sample and 3% in the pure IRAC sample. 
Both ratios are much lower than similar ratios («20%) 
in low-mass star forming re gions such as Taurus and 
Ophiuchus (jChen et al.lll995h . This might indicate that 
the star formation activity in the immediate vicinity of 
NGC 2244 is coming to an end. On the other hand the 
position of the class I sources implies that star formation 
is still ongoing outside the evacuated region . This is in 
good agreement with the results of (|2005l ) who found 
that sequential formation is taking place in the ambient 
molecular cloud. 

5. CONCLUSIONS 

We have conducted an infrared imaging survey at 3.6, 
4.5, 5.8, 8.0 and 24^tm of about 0.4 square degrees cen- 
tered on the main OB core of NGC 2244. The limiting- 
magnitudes of our photometry are 17.0, 16.4, 14.2, 13.6 



and 10.7 at 3.6, 4.5, 5.8, 8.0 and 24/im respectively. From 
the analysis of the images and the photometry we con- 
clude the following. 

• There is a strong 24 /im feature in the middle of 
the main OB core of NGC 2244. It has bowshock- 
like geometry probably shaped by the stellar wind 
of the nearby O and B stars interacting with the 
cool dust in the neighborhood. 

• We identify 337 class II sources and 25 class I 
sources on the IRAC two color diagrams and 213 
class II and 20 class I sources using the IRAC 
datasets combined with 24 /im observations. 

• In both samples the class I/class II ratio is much 
lower than similar ratios in low-mass star-forming 
regions, implying that the star formation activity 
in the immediate vicinity of NGC2244 is coming to 
an end. 

• The center of concentration of the class II sources 
shows very good agreement with the c oordin ate of 
the center of the cluster derived byO (|2005l ) using 
star-counts from the 2MASS database. 

• Using the IRAC (and where available MIPS) data 
we find that roughly 44.5% of cluster members have 
disks, which is somewhat smaller but in reasonably 
good agree ment with the result s of lHaisch et al.l 
(|2001f ) and (jMamaiek et al.ll2004D for clusters with 
similar age. 

• We found 16 class II sources in the immediate 
vicinity of O stars. We cannot completely explain 
this by line-of-sight coincidence (the most proba- 
ble number of chance alignments is 9). However, 
it is reasonable to assume that these sources have 
disks that have not been close long enough to be de- 
stroyed by photoevaporation. They might be sim - 
ilar to the sources reported by iBalog et alj (J2006) 
but unresolved because the tail is not on the plane 
of the sky. About 2/3 of these object (10-11 in our 
case) remain unresolved at 24 /im due to the large 
beamsize of Spitzer. 

• We calculated the distance of the class II sources 
from the nearest O star and from the corrected dis- 
tribution we conclude that the fraction of stars with 
disks does not correlate with the distance outside 
0.5 pc from the O stars. However there is a deficit 
of disks within the 0.5pc radius region. 

• For disks beyond 0.5pc from the O stars, we showed 
that the measured disk fraction favors the steeper 
slopes (a > 1) in the photoevaporation rate vs 
O star distance rel ation; that is slop e s like the 
a = 0.767 value of iRichling fc Yorkel (jl998l) are 
unlikely in the case of NGC 2244. This together 
with the small difference (< 10%) between the disk 
fraction of clusters with and without O stars and 
the fact that the disk fraction is constant outside a 
certain radius from the O stars, indicates that the 
effect of high mass stars is significant only in their 
immediate vicinity (d < 0.5 pc). 
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• Allowing for projection, it is likely that all of the 
identified class I sources are located around the 
rims of the surrounding nebula. This indicates 
that they arise from self-propagating star forma- 
tion triggered by the expanding shell driven by the 
high mass stars. It also implies that the class I 
objects are extremely young. 
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